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Molecular dynamics simulations with separate thermostats for rotational and translational motions were used 
to study the effects of these degrees of freedom on the structure of water at a fixed density. To describe water 
molecules, we used the SPC/E model. The results indicate that an increase of the rotational temperature, Tr, 
causes a significant breaking of the hydrogen bonds. This is not the case, at least not to such an extent, when the 
translational temperature, Tj, is raised. The probability of finding an empty spherical cavity (no water molecule 
present) of a given size, strongly decreases with an increase of Tr, but this only marginally affects the free 
energy of the hydrophobe insertion. The excess internal energy increases proportionally with an increase of Tr, 
while an increase of Tj yields a much smaller effect at high temperatures. The diffusion coefficient of water 
exhibits a non-monotonous behaviour with an increase of the rotational temperature. 
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1. Introduction 

We dedicate this contribution to J.A. Barker and D. Henderson m. authors of the important paper 
What is “liquid”? Understanding the states of matter, which was published in Reviews of Modern Physics, 
almost exactly 50 years ago. This excellent paper, which served as a textbook to many of us, is still a ref¬ 
erence article for every person who wishes to study the theory of liquids. Our contribution here touches 
upon one particular, but extremely important liquid — water, from the standpoint of the individual de¬ 
grees of freedom. 

Complex molecules have, besides the translational, also rotational and vibrational degrees of free¬ 
dom. They all make their own contribution to thermodynamic functions, which depend on temperature 
of the system. Temperature is accordingly considered to be the most important thermodynamic param¬ 
eter. An interesting question is how the fundamental degrees of freedom contribute to the potential of 
mean force between water molecules, when their characteristic temperatures vary independently. 

A biologically important situation arises when water molecules are irradiated with microwaves of 
appropriate frequency to excite their rotational motion, which affects the solvation and interaction be¬ 
tween the solute molecules (2}il. The interaction of microwaves with water has been previously studied 
using the non-equilibrium molecular dynamics simulation O, where the electric field was modelled ex¬ 
plicitly. It has been shown that microwaves induce heating by excitation of rotational motion in water 
molecules. The excess rotational kinetic energy is then transferred to translational degrees of freedom. 
Following previous studies |2}0l, the direct interaction of microwaves with a solvent is replaced by a 
model system in which the rotational degree of freedom has a different temperature than the transla¬ 
tional motion. Such a non-equilibrium steady-state picture, where different degrees of freedom can have 
different temperatures, has already been used to explain experimental phenomena, such as increased 
chemical reaction rates and altered polar solvent boiling points ElO. 

In the preceding contribution |4], we studied the effects in hydration of cations, anions and hy¬ 
drophobes caused by independent variations of the rotational and translational temperatures. The main 
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conclusions were: (i) an increase of the rotational temperature affects the hydration of cations in an op¬ 
posite way to anions, (ii) an increase of the translational temperature always decreases the height of the 
first peak in the solute-water radial distribution function; (iii) an increase of the rotational temperature 
yields an increase in the first peak in the solute-water radial distribution function for hydrophobes and 
cations; while, (iv) in contrast to this, the solvation peak decreases around ions with a sufficiently large 
negative charge. 

The focus of the preceding work was, therefore, on solvation of various solutes and not much atten¬ 
tion was paid to the solvent itself. With this respect, the present work complements the previous one; 
here, we are primarily interested in the structure of the water under such non-equUibrium conditions. In 
order to explore the properties of the model water, we used molecular dynamics (MD) simulations, where 
the rotational and translational motions were coupled to separate thermostats. We used a rigid model of 
water (SPC/E) Q; the effect of vibrational temperature variations was found to be smaU (S). The focus of 
the present study is on the structure of water as reflected in various site-site distribution functions. The 
probability of water molecule to form a certain number of hydrogen bonds is calculated as a function of 
rotational and translational temperatures. The effects of the separate Tr and Tj variations on the excess 
internal energy of the system and on the dynamics of water molecules is examined. 


2. Model and simulations 

Molecular dynamics simulations were performed at a constant number of particles and volume (AT, V) 
with 256 SPC/E water molecules in the simulation box. As before (4), the number density of water mole¬ 
cules was set to 1.0 g/cm^, and the time step used was 1.0 fs. Each MD simulation had at least 100 ps long 
MD equilibration run, while the statistics was collected during the 4900 ps long production run. 

The translational motion of a single water molecule i was described by its center-of-mass position r,-, 
velocity v,- and the force f; acting on it. Similarly, the rotational motion of a single water molecule i was 
characterized by its orientation (using quaternion representation) q,, angular velocity oji and torque Tj. 
While the acceleration of the center-of-mass follows immediately from the force acting on it, the angular 
acceleration must be computed from Euler’s equations of motion for rigid bodies. We employed the stan¬ 
dard Verlet algorithm to integrate the equations of motion (^. To hold Tr and Tj fixed at their respective 
values, we employed a simple velocity re-scaling scheme. 


3. Results and discussion 

3.1. Radial distribution functions 

We start the presentation of numerical results with the oxygen-oxygen, gooi and oxygen-hydrogen 
distribution functions, goH- Eirst, in figure]^ we present the goo distributions for various rotational (Tr) 
and translational (Tx) temperatures of the system. 

In both cases, an increase of the temperature decreases the height of the first peak of the oxygen- 
oxygen distribution function. However, there is a difference regarding how the position of the peak 
responds to the temperature variations. An increase of Tr causes the first peak to move toward larger 
distances. In contrast to this, the position of the first peak in goo does not change upon an increase of 
the translational temperature, Tx. In addition, the peak becomes somewhat broader; the shape seems to 
reflect the fact that molecules with larger kinetic energy can easier “penetrate” into each other. 

Liquid water is at ambient conditions characterized by the position of the second peak in this figure 
(figure [^. The peak is located at 4.5 A rather than at ~ 6 A as it is at the absence of hydrogen 

bonds. The position of this peak is a consequence of the hydrogen-bonding of water molecules. The shift 
of the second peak toward ~ 6 A suggests that some hydrogen bonds break upon an increase of rotational 
temperature. The effects of temperature variations on the hydrogen bonds in the model water are most 
clearly shown in figure]^ Here, the oxygen-hydrogen distribution, gon. is presented as a function of the 
respective temperature. The position of the first peak of gon distribution is characteristic of a hydrogen 
bond. We see that the height of this peak is very sensitive to Tr and much less to the Tj variations. In both 
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Figure 1. (Color online) Oxygen-oxygen radial distribution functions (goo) for different Tr (left-hand) 
and Tj (right-hand) values. Different lines correspond to different values of the relevant temperature: 
300 K (solid, red), 500 K (dashed, green), 700 K (dash-dotted blue), and 900 K (dotted, magenta). 




Figure 2. (Color online) Oxygen-hydrogen radial distribution functions (goe) for different Tr (left-hand) 
and Tj (right-right) values. The color code as for figure[^ 


cases, however, the peak decreases in magnitude. From this distribution, the coordination number, uh. is 
calculated (for the definition see caption in figure]^. This quantity, as shown in figurej^ decreases almost 
linearly with Tr. On the other hand, tin is much less sensitive to the Tj variations at higher temperatures. 



T[K] 


Figure 3. (Color online) Coordination number, wh. of hydrogen atoms of the neighbouring molecules 
around oxygen of the central molecule (and vice versa) as a function of Tr (solid, red) and (dashed, 
green). «h is obtained as the volume integral of gou up to the first minimum (at 2.4 A). 
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Figure 4. (Color online) Probability for a water molecule to form Whb hydrogen bonds, PIWhb). at differ¬ 
ent Tr (left-hand panel) and Tj (right-hand panel). We count two water molecules as hydrogen-bonded if 
(i) the angle between the OH arm of one molecule and the 0-0 vector is less than 30°, and (ii) if the 0-0 
distance is less than 3.3 A (position of the first minimum in goo)- The color code as for figurej^ 


In figure]^ we show the probability PCAThb) that a water ntolecule forms A/hb hydrogen bonds. Up 
to 500 K, the probability distribution, P(A(hb). is practically the same for the Tr and Ti dependence. 
Above rR=500 K, the peak of this distribution is gradually shifted toward smaller A/hb values and there 
is an appreciable probability for only one hydrogen bond. The situation is different for an increase of Tj 
where the PCA/hb) distribution does not change much upon an increase of the translational temperature 
above rx=SOO K. 

All the evidence gathered so far suggests that hydrogen-bonding is more sensitive to Tr, than to Tj 
variations. The reason for this has been discussed in reference (3) (page 024108-6). The formation of a 
hydrogen bond requires two water molecules to meet at an appropriate distance and a correct orienta¬ 
tion. It is easy to understand that every increase in Tr (faster rotation) whl make favorable orientations 
less probable. At conditions of high Tr, the model molecules will behave as independent rotators with 
the hydrogen bond interaction energy approaching zero. 

The mechanism causing two molecules of water to lose the correlation upon the Tj increase is differ¬ 
ent. In the latter case, the molecules should separate for a certain distance from each other, which is not 
that easy because each molecule is caged by the others. By increasing Tj, having in mind that we perform 
simulations at a constant water density, we do not approach the ideal gas behaviour, even at high tem¬ 
peratures (rT=900 K) the system is under the effect of strong repulsive interaction. An initial increase of 
Tf (up to Tj - 500 K) indicates that some hydrogen bonds are broken and the radial order is disrupted to 
a certain extent. 

A more detailed information on the water structure can be obtained if the radial distribution function, 
goo. is divided into the separate contributions of the neighboring molecules coordinated on the chosen 
molecule. In our notation, gj is the contribution of the i-th neighbor to the total pair distribution function. 
These results are shown in figure where the left-hand panels belong to different rotational and the 
right-hand panels belong to different translational temperatures. In case of Tr increase, we observe the 
shifts of the peaks toward larger distances. This applies to all four nearest neighbors. At this density water 
molecules are closely packed; while aww for this (SPC/E) model is 3.169 A, the position of the first peak 
is actually at around 2.7 A. Rising the Tr causes for this distance to increase, and the molecules distribute 
at larger average distances. In contrast to this, an increase of the Tj does not significantly change the 
position of the first peak. This is yet another consequence of the fact that faster rotation (increase in Tr) 
disrupts the hydrogen bond network more efficiently than an increase in Tj. 

An interesting question, related to the changes of peak positions of various g,’s upon Tr increase, has 
been posed by an anonymous reviewer. Figure [^indicates that the position of the first peak, gi, is less 
affected than the positions of more distant peaks, such as ga and g 4 . Why? Considering that ga measures 
the contribution of the third neighbour to the oxygen-oxygen distribution function, we may conclude 
that correlation between the central molecule and the third sheU molecules is significantly suppressed. 
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Figure 5. (Color online) Oxygen-oxygen radial distribution function for the first four neighbours of the 
water molecule (from top to the bottom gi to g^) at different Tr (left-hand panels) and Tj (right-hand 
panels). The color code as for figurej^ 


In other words, the effect of the central molecule does not propagate that far in the bulk as for T^-Tj- 
300 K. This is clearly shown in figurej^ the structure of the solution as measured by oxygen-oxygen pair 
distribution function, goo. becomes more uniform upon an increase of either Tr or Tj. Furthermore, 
oscillations are not just smaller at elevated temperature, they are also “out of phase”. 
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3.2. Probability of observing an empty cavity of a size of water 

The simulations presented in this paper are performed for a constant number of particles in a given 
volume. An important information on a liquid system is associated with the density fluctuations. In this 
subsection we studied the probability distribution P{N) for observing exactly N water molecules (ac¬ 
tually their centers of gravity) in a spherical volume having diameter 3 A. This probability distribution 
is approximately Gaussian (9] for smaU (molecular size) volumes. The width of the distribution reflects 
compressibility of the liquid: a wider distribution means a somewhat higher compressibility. Our sim¬ 
ulations indicate that P{N] curves change only marginally upon an increase of Tr and Ti- A notable 
exception is the P(0) value — i.e., the probability of observing an empty spherical cavity. Accordingly, 
only the dependencies of P(0) on Tr and Tj are shown in figureAn increase of Tr from 300 to 900 K 
causes a decrease for about one order of magnitude for P(0) (note the natural logarithm plotted on y-axis 
in figurej^. The density fluctuations are suppressed due to a decreased number of hydrogen bonds. On 
the other hand, the variations of Ti have little effect on P(0). 



T[K] 


Figure 6. (Color online) Logarithm of probability P(0) of observing an empty spherical cavity of diameter 
3 A in the model water as a function of Tr (solid, red) and (dashed, green). 

Complementary to these calculations, we also performed the thermodynamic integration to simulate 
the free energy of a hydrophobe insertion, APms, into the model water (see figure]^. Hydrophobe was 
represented by the Lennard-Jones particle with the well depth equal to that of the model water molecule 
and ass = 3 A. Lorentz-Berthelot mixing rule was used for cross interaction parameters. Under conditions 
where the hydrophobe-water coupling is weak, we expect for the APins to have qualitatively the same 
features as the same quantity for the hard-sphere solute. 

While the P(0) values decrease with an increase of Pr, the insertion free energy is only slightly af¬ 
fected under such conditions. On the other hand, Ti, which has virtually no effect on the P{N) distribu¬ 
tion (see figure]^, has a relatively large effect on the APjns. In both cases, the free energy for the particle 
insertion, APi^s, increases; it becomes increasingly more difficult to add a hydrophobe to the system. 

Here, we offer a qualitative explanation for this result. The insertion free energy of the hard-sphere 
solute is related to the probability of observing an empty spherical cavity, P(0), in solvent @ (as usual 
is the Boltzmann constant): 


APins = -A;BrinP(0). (3.1) 

As the solute-water interaction depends only on the positions of water molecules (and not on their ori¬ 
entations) and since equation (13 applies to any solvent (not necessarily possessing rotational degrees 
of freedom), we might expect that temperature in equation (13 refers to the translational temperature. 
While the P(0) values vary substantially with Pr, the InP(O) only increases for about 20 % in magnitude, 
which is close to the relative change of APins with Pr. Much stronger increase of AP^s is observed upon 
the rise of Tj and it may be attributed to an increase of the pre-factor T. 
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Figure 7. (Color online) Free energy of the hydrophobe insertion as a function of Tr (soUd, red) or Tj 
(dashed, green). 


3.3. Excess internal energy 

It is of interest to examine the variations of the excess internal energy (Eex) per mol of water molecules 
as functions of Tr or Tj (figurej^. We can see that elevation of Tr causes an almost linear increase of itex. 
while the same quantity changes less and less if Tj is increased. The reason for such a behaviour lies, as 
explained before, in different effects of the two degrees of freedom to hydrogen bonding. These graphs 
also provide some insights into the Tr and Tj effects on the excess heat capacities (at constant N and V), 
C“. and . As we see, the rotational contribution (derivative of the upper curve) is temperature inde¬ 
pendent, while for the translational contribution (derivative of the lower curve) the excess heat capacity 
goes to zero for high Tj values. 

3.4. Dynamics of model water molecule upon an increase of Tr and Tj 

How the different degrees of freedom affect the dynamics of water molecules? We can see the answer 
in figure[9] showing the effect of Tr and Tr on the (self) diffusion coefficient of a model water molecule. 

At first, an increase of Tr causes a faster diffusion, which can be attributed to the breaking of hy¬ 
drogen bonds. As the water molecules are more weakly bound, they tend to diffuse faster. However, an 
increase of Tr also means a decrease of P(0). Therefore, a random walk should proceed in smaller steps, 
causing the diffusion to slow down. This effect seems to dominate at the highest Tr studied here. On the 
other hand, increasing Tx. as expected, causes a monotonous increase of the diffusion coefficient in this 
temperature interval. 



Figure 8. (Color online) The excess internal energy per mol (i?ex) as a function of Tr (soUd, red) and Tj 
(dashed, green). 
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T[K] 

Figure 9. (Color online) Diffusion coefficient as a function of Tr (solid, red) or Tj (dashed, green). Note 
the logarithmic scale on y-axis. 


4. Conclusions 

It is well known that microwaves induce heating by excitation of rotational motion in water molecules 
and that the rotational kinetic energy is then transferred to the translational degrees of freedom. The 
microwave irradiation is important in biology: there are reports that microwaves enhance the folding 
and unfolding kinetics of globular proteins I10M2I . as well as protein aggregation (13). 

The effect of microwaves on the properties of liquids and solutions can be studied by non-equilibrium 
molecular dynamics simulations O. In such studies, the electric field is modelled explicitly. Following the 
ideas, proposed in references (2] [3), we perform calculations in which the rotational degree of freedom 
has the temperature different from the translational one. In a previous paper (4|, we examined the effect 
of such non-equilibrium conditions on the solvation of simple solutes. In the present contribution, we 
studied the properties of a model water under non-equilibrium conditions where the translational and 
rotational temperatures vary independently. The main conclusions of this work are: (i) an increase of the 
rotational temperature causes the first peak of goo to move toward larger distances, while its position 
remains unchanged upon Tj increase, (ii) The height of the first peak in the hydrogen-oxygen distribution 
function is much more sensitive on the Tr than on the Tj variations, (hi) At high rotational temperatures, 
the number of hydrogen bonds may drop as low as to one. (iv) Free energy of the hydrophobe insertion 
does only marginally depend on the rotational temperature, (v) Water diffusion coefficient measured via 
the mean square displacement is, in contrast to the Tj dependence, a non-monotonous function of the 
rotational temperature. These and other effects observed upon varying the rotational temperature may 
significantly affect the properties of water as solvent, which has to a certain degree been examined in 
reference o The work exploring such effects on polymer and poly electrolyte conformations in aqueous 
solution is currently underway — preliminary results clearly point toward a high probability of polymer 
coUapse in cases of rotational heating. 

In the end, we shall conclude that, despite important developments following reference (1), there is 
stUl ample room for improvements in theory and simulation of liquids, of which, water is by far most 
important. 
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Bn/ii/iB TpaHC/iJii^iiiHi/ix ra o6epTOBi/ix CTyneHiB Bi/ibHOCTi Ha 
B/iacTi/iBOCTi MOAe/ii boai/i 

T. Moxopii, B. rpi6ap-/li, B. B.iiaxi 

OaKyabteT xImIT I xiMlHHOiTexHoaorii, VHlBepcnTeT /IroSanHH, AutKepneea 5, SI-1000 /IK)6anHa, CaoBeHin 

flaa BUBneHHa BnanBy rpaHcaaqlMHnx ta oOeproBnx cryneHlB BiabHOcri Ha crpyKtypy BOfln 3 cJlxcoBaHOK) ry- 
CTi/iHOK) BHKopMCtaHO MBTOfl MoaeKyaapHOi flHHaMlKH 3 flBOMa OKpewi/iMH TepMOCTaraMi/i flaa o6epTOBnx ta 
rpaHcaaqlHHi/ix pyxlB. flaa onncy Moaexya BOfln 6yao BHKopucraHO MOfleab SPC/E. PeayabTan/i noKaayiOTb, tqo 
36iabUjeHHa o6epTOBOi teMneparypH, Tr, npn3BOAMTb flo 3HaHHOro po3pnBy BOflHeBi/ix 3B'a3KiB. IJiaKOM 3bo- 
poTHlM e BHnaflOK, npHHaitMHi He flo Taxoi Mipit, Koan aSiabiuycTbca rpaHcaaqlitHa reMneparypa, Tj. klMOBip- 
HiCTb 3HaHTi/i ccJepHHHy nopoxHHHy (BiflcytHa Moaexyaa BOfli/i) aaflanoro po3Mipy ananHO 3MeHUjyeTbca npi/i 
36iabUjeHHi Tr, aae qe Mae amue HeanaHHHH BnawB na BiabHy eneprlK) riflpoc|)o6HoT BcraBKia. HaflawujKOBa 
BHyipituHa enepria 36iabUjyeTbca nponopqlHHO flo 36iabUjeHHa T^, b tow aac, ax Tj flae anaano MeHiui/iH 
ecJeKT npH bmcokhx reMneparypax. KoecJiqieHT flnc|)y3iT BOflia Mae HeMOHOtOHHy noBe^lHKy npi/i 36iabUjeHHi 
o6epTOBoTTeMnepaTypi/i. 

K/iioaoBi cnoBa: crpyKiypa BOfl,i/i, CTyneni Bi/ibHOcri, MoneKynupna fli/inaMiKa 
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